Coherent current transport in wide ballistic Josephson junctions 
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We present an experimental and theoretical investigation of coherent current transport in wide 
ballistic superconductor-two dimensional electron gas-superconductor junctions. It is found exper- 
imentally that upon increasing the junction length, the subharmonic gap structure in the current- 
voltage characteristics is shifted to lower voltages, and the excess current at voltages much larger 
than the superconducting gap decreases. Applying a theory of coherent multiple Andreev reflection, 
we show that these observations can be explained in terms of transport through Andreev resonances. 
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Ballistic superconducting junctions are of great inter- 
est for studying fundamental properties of coherent elec- 
tronic transport. The concept of Josephson effect in 
superconductor-normal metal-superconductor junctions 
was originally formulated for ballistic junctions 1]. How- 
ever, until recently, experiments could only be performed 
on metallic diffusive junctions, which has limited the 
range of phenomena observed. Experimental realiza- 
tion of the ballistic regime became possible when a high- 
mobility two-dimensional electron gas (2DEG) was em- 
ployed to connect superconducting (S) electrodes 0- 
The recent experimental interest [j| has been focused on 
junctions with InAs 2DEGs which form highly transpar- 
ent 2DEG-S interfaces with large probabilityof Andreev 
reflection. Moreover, electrostatic gating 0,0 has made 
it possible to control junction parameters. 

The dc Josephson effect in ballistic S-2DEG-S junc- 
tions has been extensively studied both experimentally 
0. El & El and theoretically @] . Less investigated are 
the transport properties of the junctions in the presence 
of a voltage bias between the superconductors, the regime 
of multiple Andreev reflection (MAR) 8J. In particular, 
a conclusive experimental picture of the subharmonic gap 
structure (SGS) in the current-voltage characteristics is 
lacking. Moreover, the main tool for analyzing the ex- 
perimentally observed current- voltage characteristics has 
been the theory of Octavio et al. |3| (OTBK). However, 
this is not appropriate for junctions which show a dc 
Josephson effect, i.e. which are in the coherent transport 
regime, since the OTBK theory is restricted to incoherent 
MAR transport. 

The theory for coherent MAR transport in junctions 
much shorter than the superconducting coherence length 
flf)j | has with great accuracy explained experiments with 
atomic point contacts both the details of the SGS at 



eV = 2 A jn and the magnitude of the excess current at 
voltages much larger than the superconducting gap. This 
theory however cannot be directly applied to the exper- 
imentally studied junctions 0- lU 0- E3|- generally of 
lengths L of the order of the superconducting coherence 
length £ - 

In this letter we report a systematic experimental in- 
vestigation of coherent current transport in ballistic Nb- 
InAs-Nb juctions with 2DEG-S interfaces of high trans- 
parency. For every junction we investigated the three key 
transport characteristics: SGS, excess current, and criti- 
cal Josephson current. All investigated junctions show a 
dc Josephson effect. Our main experimental observation 
is that the SGS is shifted down in voltage upon increas- 
ing the distance between the superconducting electrodes, 
i.e. the junction length. Moreover, the excess current is 
found to decrease with increasing junction length. 

These observations can be explained within a coherent, 
multimode MAR-theory. The shift of the SGS results 
from transport through Andreev resonances fl2l : the 
Andreev resonances, closely related to the Andreev levels 
carrying the dc Josephson current, are shifted to lower en- 
ergies upon increasing the junction length, consequently 
shifting the SGS down in voltage. Similarly, the excess 
current decreases with increasing junction length for the 
experimentally relevant junction lengths, L < £n. 

A schematic picture of the junction is shown in Fig. 
H Two niobium electrodes are deposited directly on the 
2DEG which forms spontaneously at the surface of single 
crystal InAs I n the following step, the junction is 

covered by a Si02 layer that provides insulation to the 
Al gate electrode grown on top. On three different chips 
A,B,C, junctions of five different interelectrode distances 
L = 105(A), 115(B), 145(B), 160(A) and 200(C) nm have 
been fabricated. All junctions have a width W = 50 /im. 
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The sheet carrier density n a , the effective mass m* and 
the mobility \x are controlled with the top gate to give 
the largest 2DEG-conductance. Typical parameters for 
the unprocessed InAs are n s = 2 x 10 16 m~ 2 and m* = 
0.033m o , where m c is the free electron mass. This gives a 
Fermi wave length Af = 18 nm and vf = 1.2 x 10 6 ms _1 . 
The mobility was not measured on the actual samples, 
but is typically of order fi = 10000 cm 2 (Vs) -1 . The 
measured critical temperature of Nb T c = 9.0 K corre- 
sponds to a superconducting gap A = 1.4 meV close to 
the bulk value. The superconducting coherence length is 
thus £o = TiVf/A — 600 nm. The normal resistance J?n 




FIG. 1: Left: Average normal transparency D = R^/Rsh as a 
function of n s for samples of lengths L — 105 nm and 160 nm. 
Right: I c Rn product as a function of T for the same samples. 
The thin lines are the theory, with I C R^(T = 0) scaled for 
best fit by 1.65 (L = 105 nm) and 2.4 (L = 160 nm) and 
Tivf/Ea = 0.54£o (see text). Inset: A schematic cut through 
the junction, displaying the kink in the 2DEG. 

was measured at a voltage eV 3> A. The ratio between 
the Sharvin resistance i?sh = h/ (2e 2 )[\F /2W] and the 
normal resistance, i.e. the average transmission proba- 
bility per conduction mode, is D ~ 0.7 for n s > 2 x 10 16 
m -2 (see Fig. QJ. The interfaces between the 2DEG and 
the metal electrodes are thus highly transparent. At- 
tributing the small difference in resistance between the 
105 and 160 nm junction to impurity scattering gives a 
mean free path I 300 nm. This demonstrates that 
the transport is in the ballistic regime (L < I). All 
junctions show a dc Josephson current, with I c i?N prod- 
ucts 0.6,0.55,0.52,0.45 at 1.8 K monotonically decreas- 
ing with L from L = 105 to 200 nm. The temperature 
dependence of the dc Josephson current is shown in Fig. 
^for two junctions L — 105 and 160 nm. 

In Fig. [2] the experimentally measured differential re- 
sistance for the different junctions is depicted. For all 
junctions, the resistance shows SGS at voltages eV < 2A. 
As an overall tendency, the SGS is shifted to lower volt- 
ages when increasing junction length. For junctions on 
the same chip this holds strictly. All SGS-features are 
symmetric around V = and scale roughly with temper- 
ature like the gap A(T), up to T c . 

A decrease in the excess current as a function of length 
L is also observed in the experiment. The 7 exc i?N prod- 
uct for the different junctions is shown in Fig. 21 below, 



displaying an overall decrease with increasing L. This 
holds strictly, just as the shift of the SGS, for junctions 
on the same chip. The systematic shift of the SGS and 
the decrease of the excess current with junction length 
are the main experimental observations in this work. 

Based on these findings, we model the junction as a 
ballistic, two-dimensional normal conductor of length L 
and width W connected to two bulk superconducting 
electrodes via highly transparent interfaces. The nor- 
mal reflection in the junction is assumed to be concen- 
trated to the kink in the 2DEG formed during the junc- 
tion processing, see Fig. ^ For simplicity, the normal 
reflection is taken to be specular, reducing the problem 
to summing over independent modes. The width of the 
junction W ^> Af, giving a number of transport modes 
2W/Xf = 5600. No sources of decoherence, e.g. inelastic 
scattering [lfij ]. are taken into account in the model. 

We perform our theoretical analysis within the scat- 
tering approach to the Bogoliubov-de Gennes equation 
[lif . This allows us to express the dc Josephson cur- 
rent 1171 as well as the current in voltage biased junc- 
tions 12] in terms of the transmission and reflection 
amplitudes t m (E) and r m (E) of each transport mode 
m and the mode independent Andreev reflection prob- 
ability a(E) at the 2DEG-S interface. The scattering 
amplitudes are given by r m = iy/R(l + q m )/(l + R<lm) 
and t m = q m (l — R)/(l + Rq m ) where q m = exp[ik m L]; 
h m (E) = [2m*/h 2 {E F - E m + £)] 1/2 , where E F is the 
Fermi energy and E m is the transverse mode energy. 
The reflection probability of each interface R, taken to 
be mode independent, is given from the averaged total 
transmission probability D as R = (1 — D)/(l + D), giv- 
ing R = 0.16 for D = 0.72 (D for the 105 nm junction at 
n s = 2 x 10 16 m~ 2 ). 

The Andreev reflection amplitude at the 2DEG-S in- 
terface is given by a{E) — exp[— iarccos(£?/A)]. It is 
however not possible to apriori conclude that quasipar- 
ticles are Andreev reflected directly at the 2DEG-S in- 
terface near the kink, they might spend a certain time 
h/Ej^ underneath the superconductor [l^. This can be 
incorporated into the model in the simplest possible way 
by multiplying a(E) with a phase factor exp(iE/ E&), 
where E& has to be determined by comparison to the 
experimental data. 0] 

The total dc-current can be written as a sum of the cur- 
rents for the individual modes [2(J, i.e. I = ^2 m I m - In 
the general case, it is not possible to analytically calculate 
the current for arbitrary voltage, but one has to resort 
to numerics. It is important to note that most transport 
modes m have effective lengths L/^l — (m/N) 2 close to 
the physical length L. Thus, the sharp features in the 
SGS predicted in single mode junctions of finite length 
|l3| are not washed out but are merely broadened by the 
summation over the modes. 

The qualitative behavior of the SGS as a function of 
junction length is clear from the numerics in Fig. |3J One 
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FIG. 2: Left: Measured differential resistances dV/dl (upper) 
and (numerically differentiated) d 2 I/dV 2 (lower) as a function 
of V for T = 1.8 K. The junction lengths are from bottom 
to top, L = 105, 115, 145, 160 and 200 nm. The traces have 
been successively shifted for clarity. Right: Positions of the 
minima Vmin of the d 2 I/dV 2 (lower left) for the SGS n = 1 
(<), n = 2 (•), n = 3 (o), n = 4 (l>, only two minima clearly 
visible) and n = 5 (o). Dashed lines are given by Eq. Q, 
with an additional length Hvf/Ea = 0.6£o- 



sees how the SGS are shifted down to lower voltages when 
increasing the junction length L. The nature of this ef- 
fect is the shift of broadened MAR resonances 0, 0] . 
The MAR resonances are closely related to the Andreev 
levels carrying the dc Josephson current. The Andreev 
levels are shifted towards lower energies upon increas- 
ing the junction length, leading to a corresponding shift 
in the SGS. Additional numerics shows that the general 
tendency of a shift of the SGS towards lower voltages 
is independent of the reflection probability R, while the 
details of the SGS depend on R. For perfectly transpar- 
ent 2DEG-S interfaces, however, the SGS are absent for 
fcT«A| for any junction length [l3| . 

Quantitatively, we have performed a detailed numeri- 
cal investigation of the SGS (n = 1 — 5) in the experimen- 
tally relevant limit of small normal reflectivity, i? < 0.16 
and small effective junction length, L c h < £o- The ef- 
fective length L c g — L + Tiv-p / E\ is the sum of the geo- 
metric length of the junction L and the distance Tivf/Ea 
the quasiparticles propagate under the superconductors 
before they are Andreev reflected. The analysis of the 
numerical results shows that close to V = 2 A /en the sec- 
ond derivative [2l| d 2 I/dV 2 has minima which depend on 
junction length as 



V m 



2A 



1 - 0.4- 



1,2,. 



(1) 



for L c ff < 0.4£o for all n (for n > 1, it holds for even 
longer L e ff) ■ This length dependence is in strongcontrast 
to the result of the theory of incoherent MAR |9j , which 
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FIG. 3: Left: Plot of the numerically calculated dl/dV as 
function of eV/A and L/£ for R = 0.16 and kT < A; 
light areas correspond to peaks, dark to dips. (Tlv-f/Ea = 0). 
Right: A contour plot of the third derivative dl 3 jdV 6 around 
zero, giving the extremal points of dl 2 /dV 2 in the intersec- 
tions of the dark and white areas. The minima in Eq. Q arc 
indicated by dashed lines. 



predicts dips in the second derivative at voltages V m - U1 = 
2 A/en independent of junction length p^ . 

Comparing theory and experiment, we see that the co- 
herent MAR theory (Fig. [3J is able to explain the overall 
shift in the experimentally observed SGS towards lower 
voltage for increasing junction length L (see Fig. 
However, to fit the measured positions of the dips in the 
second derivative with Eq. (JTJ, one needs to include an 
additional length fiv-p/E^. The best fit, including the 
SGS to fifth order, is obtained for Hvf/Ea = 0.54£ for 
both junctions on chip A, 0.7£o for both junctions on 
chip B and 0.7£o for chip C. We note that the values for 
the effective lengths L c s are partly outside the range of 
validity of Eq. (1), which might introduce some error in 
the fitted values Hvf/Ea- The fact that the junctions 
on the same chip give the same Ea suggests that E\ is 
sensitive to the sample preparation. 

Although the dip and peak positions of the SGS can 
be well reproduced by the theory, the theoretically cal- 
culated amplitudes (not shown) of the SGS-features are 
much larger than the ones observed in the experiment. 
One reason might be residual inelastic scattering or de- 
phasing, not taken into account in the model. Another 
reason might be that the theory assumes a voltage bias, 
while the measurements are carried out with current bias, 
giving rise to voltage fluctuations which smear the SGS. 

The experimentally measured excess current I CKC = 
I — V/Rn at voltages V ~> A/e is presented in the 
left panel in Fig. As seen from numerics in the 

right panel of Fig. H for both E A = 1.4 A and 1.8 A, 
the excess current decreases monotonically with junction 
length L for all lengths L < £o, in qualitative agreement 
with the experiment. However, this decay is slower than 
the measured one. Moreover, the theoretical excess cur- 
rent exceeds the experimental one. This discrepancy can 
be attributed to residual inelastic scattering, similar to 
the suppression of the amplitude of the SGS mentioned 



4 




FIG. 4: Left: Excess current at T = 1.8 K for the five exper- 
imentally investigated junctions. Right: Numerically calcu- 
lated excess current as a function of length L, for R = 0.16 
and k B T < A. 

above. The excess current is carried by quasiparticles 
which traverse the junction one and two times. The 
single-particle processes give a negative contribution to 
the excess current while the two-particle processes yield 
a positive contribution. Since quasiparticles carrying the 
two-particle current spend longer time in the junction, 
inelastic scattering leads to a stronger suppression of the 
two-particle current compared to the one-particle cur- 
rent. This results in an excess current with a stronger 
length-dependent decay than predicted by the coherent 
theory. 

We also compare the experimentally measured dc 
Josephson currents with the result of the model. The 
expression for Josephson current in terms of the normal 
scattering amplitudes and the Andreev reflection ampli- 
tudes is given in Ref. 0. Using the value of Ea de- 
duced from the SGS fitting, we see in Fig. ^ that the 
temperature dependence of the experimentally measured 
critical current is well reproduced by the theory. How- 
ever, the theoretical overall magnitude is 1.65 times (for 
L = 105 nm) resp. 2.4 times (for L — 160 nm) larger 
than the experimental value. This discrepancy indicates, 
in line with the discussion for the SGS-amplitude and the 
excess current above, the presence of some decoherenece 
mechanism. Note that the discrepancy is much smaller 
than previously reported 0. 0, 0; because the longer 
effective junction length reduces the theoretical critical 
current to a value more similar to the experimental one. 
The significant part is that this effective length has been 
determined from a best fitting to the SGS. 

In conclusion, we have studied, experimentally and 
theoretically, coherent current transport in wide ballis- 
tic S-2DEG-S junctions. It is found experimentally that 
both the SGS and the excess current show a systematic 
dependence on length of the junction. We show that 
these observations can be qualitatively explained within 
a coherent theory of MAR. Our investigation also points 
towards additional mechanisms, e.g. decoherence, to be 
included in the theory to obtain a quantitative agreement 
with experiments. 
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